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ABSTRACT 


INVESTIGATION  OP  THE  UHP  MOBILE-RADIO  PROPAGATION  CHANNEL 
IN  A  TEMPERATE- ZONE  PORESTED  ENVIRONMENT 


by 

Kenneth  J.  White 


Experimental  verification  of  the  stochastic  power-spectral- 
density  theory  of  radio  propagation,  as  applied  to  a  mobile 
receiver  in  a  forested  environment,  is  presented.  Using  a 
novel  technique  to  achieve  f ractional-hertz  resolution  of 
the  diffuse  and  direct  wave  components  of  the  Doppler-spread 
received  spectrum,  the  resultant  data  taken  at  432  MHz  shows 
that  the  angle  of  arrival  of  the  diffuse  component  of  the 
incident  energy  at  the  receive  antenna  is  described  by  an 
uniform  probability  density  function  for  both  vertical,  and 
horizontal  polarizations.  Given  a  8.8  kilometer  transmission 
path  with  the  signal  penetrating  1.3  km  of  forest  just 
before  reaching  the  receiver,  the  ratio  of  coherent  (direct 
wave)  intensity  to  total  received  signal  power  is  42%  for 
horizontal  polarization  and  b%  for  vertical  polarization. 

The  excess  path  loss  above  f ree-space- loss  is  examined  and 
the  impact  of  the  results  on  existing  theories  of 
propagation  in  foliage  is  discussed. 
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INTRODUCTION 


This  decade  has  seen  a  rapid  increase  in  the  military  and 
commercial  demand  for  quality  high-capacity  mobile 
communications  in  both  voice  and  data  applications.  The  new 
technologies  in  mobile  communications  are  quickly  being 
embraced  by  the  US  Army  as  the  complex  dynamics  of  the 
modern  battlefield  and  the  rigors  of  the  army's  new  Air-Land 
Battle  doctrine  have  vastly  increased  the  demand  for 
communications .  Additionally,  the  intense  lethality  of  the 
latest  weapons  systems  dictates  that  both  the  fighting 
elements  and  their  headquarters  be  constantly  "on  the  move” 
ri].  To  add  greater  system  flexibility,  survivability,  and 
capacity,  the  army  is  supplementing  its  traditional 
terrestrial-based  systems  with  tactical  satellite  links  for 
both  fixed  and  mobile  subscribers. 

Essential  to  the  continued  development  of  quality  and 
reliable  mobile  communications  systems  for  the  battlefield 
is  the  formulation  and  verification  of  predictive 
propagation  models  for  the  tactical  environment.  The  Air 
Porce  Systems  Command  and  the  Array's  Communications- 
Electronics  Command  in  recent  years  have  been  instrumental 
in  the  military  effort  in  the  formulation  of  such  models  L2] 
-  [9]  while  Bell  Laboratories  personnel  and  others  have  been 


been  in  the  fore-front  on  the  civilian  side  [103  -  [223. 

In  an  attempt  to  assist  in  the  characterization  of  a 
mobile-radio  propagation  channel  through  a  forested 
environment  (uhich  comprises  a  dominant  portion  of  the 
potential  battlefield  in  many  of  the  areas  of  strategic 
interest  to  the  United  States),  this  paper  presents  the 
results  of  research  based  on  the  stochastic  pouer-spectral 
theory  of  propagation  [103,  [113  utilizing  a  fine-resolution 
spectrum  analysis  technique  developed  by  Campbell  [183. 

Pig.  1  shows  the  experimental  approach  taken  in  this  study. 

A  mobile  receiver  is  constrained  to  move  at  constant 
velocity  along  a  straight  stretch  of  road  uhich  has  been  cut 
through  a  homogeneous  forested  stand.  A  CH  transmitter  is 
placed  some  miles  from  the  forest  and  transmits  a  signal  in 
a  direction  orthogonal  to  the  road.  The  receiver  will 
receive  energy  via  the  "direct  wave"  that  propagates 
directly  from  the  transmitter  through  the  forest  (ray  Rl)  as 
well  as  waves  reflected  from  the  ground  (R2).  These  waves 
combine  at  the  antenna  as  the  "coherent"  field.  Energy  also 
reaches  the  receiver  from  those  waves  that  are  scattered 
from  the  trees  (rays  R3  and  R4>.  These  latter  waves 
comprise  the  "incoherent"  or  "diffuse"  field  incident  at  the 
antenna.  An  extensive  effort  has  been  made  in  determining 
the  characteristics  and  behavior  of  the  "coherent" 
intensity,  but  little  is  known  about  the  diffuse 


component.  The  nature  of  the  latter  is  of  concern  to 
theorists  and  designers  of  both  terrestrial  and  satellite 
■obile  links  £23]  and  will  provide  insight  into  the  overall 
propagation  aechanisa  of  the  aediua.  This  study  provides 
aeasureaents  of  the  coherent  and  incoherent  intensities 
resulting  froa  a  forest  propagation  path  at  432  HHz  and 
exaaines  the  pouer-spectral-density  of  the  received  signal. 
The  iapact  of  the  experiaental  results  on  the  existing 
general  theories  of  propagation  through  forests  uill  also  be 
discussed. 

FIELD  EXPRESSIONS  AND  STATISTICS 

Pig.  2  shows  an  aerial  view  of  the  aobile  receiver  traveling 
at  velocity  "Vo”  through  the  forested  stand  with  the  CH 
source  transaitting  noraal  to  the  direction  of  vehicle 
travel.  A  portion  of  the  diffuse  coaponent  scattered  froa  a 
tree  forward  of  the  vehicle  arrives  at  angle  "a”*  The 
direct  wave  will  be  received  by  the  antenna  at  the 
transaitter  frequency.  fc,  but  the  incoherent 
wave  uill  be  Doppler  shifted  to: 

f(o)  =  fc  +  < Vo/  X )cos  a  (1) 

Note  that  waves  scattered  froa  objects  forward  of  the 
vehicle  uill  be  shifted  higher  in  frequency  while  those 
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^scattered  from  behind  uill  be  shifted  down  in  frequency  with 
a  maximum  possible  shift  of  fB; 

fm  =  Vo/x  (2) 

The  bandwidth  of  the  received  signal  uill  therefore  be  2 
tines  fB.  por  example,  with  a  vehicle  velocity  of 
30  aph  and  wavelength  of  69.4  ca  (432  BH2),  the  bandwidth  of 
the  signal  at  the  receiver  uill  be  about  39  H2.  It  remains 
then  to  describe  the  statistics  of  the  signals  comprising 
this  narrow  bandwidth. 

An  electromagnetic  field  which  fluctuates  about  some  mean 
value,  such  as  the  signal  from  the  CH  transmitter  in  this 
study,  can  be  expressed  at  some  point  in  space  (eg.,  any 
point  along  the  vehicle  path  in  fig.  1)  as  the  sun  of  its 
aean  value  (coherent  contribution)  plus  its  2ero-aean 
contributions  (fluctuating  or  diffuse  component).  Using  the 
complex  envelope  representation  of  the  signal,  this  may  be 
expressed  as: 

u(r,t)  =  <u(r,t)>  +  uf(r,t)  (3) 

Considering  the  forest  to  be  a  random  medium  of  discreet 
scatterers  and  following  the  definitions  and  notation  of 
Ishiaaru  [241,  the  total  intensity,  then,  is  the  sun  of  the 


coherent  and  incoherent  field  intensities: 

IT  =  <lu(r,t)l2>  =  Ic  +  I±  (4) 

=  l<u<r,t)>l2  4-  <luf.<r,t)l2> 

where  the  brackets  indicate  an  ensenble  average  and  **tn 
indicates  a  tine  dependancy.  If  the  nediun  of  propagation 
consists  of  fixed  randon  scatterers,  then  the  propagation  is 
at  least  a  wide-sense  stationary  (and  presunably  ergodic) 
process.  Ensenble  averages  are  then  equivalent  to  tine 
averages  [25]  and  neasurenents  of  the  intensity  of  both 
conponents  of  the  signal  "u<r>”  nay  be  nade  at  a  nunber 
points  taken  sequentially  rather  than  sinultaneously.  This 
is  the  approach  used  in  this  study. 

By  treating  the  field  at  any  point  in  the  vehicle  path  as 
a  randon  variable,  the  autocorrelation  function  Ruu<x> 
between  the  field  at  any  two  points  separated  in  tine  by  "t" 
is  expressed  as: 

R^t)  =  <U<t)U#(t+T>>  (5) 

*  Lim~  /T  u(t)  u*(t+r)dt 
t-*-  ®  1  • 

Expanding  the  integrand  and  dropping  terns  which  result  in 
zero  nean  results  in: 
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Rmjtx)  =  I  <u>  I2  +  <uf<t)u£<t  +t)>  (6) 

=  Ic  +  Rf<T> 

uhere:  R^(t)  is  the  autocorrelation  of  the  fluctuating 
component. 

The  above  autocorrelations  are  expressed  in  "r”  or  the 
separation  in  tiae  between  points  being  considered.  Since 
this  separation  is  a  function  of  vehicle  velocity,  the  above 
equations  «ay  be  expressed  as  spatial  correlations  through 
the  relationship: 

distance  d  =  Wo  *  t  (7) 

Por  any  random  process  ”x<t>”,  the  autocorrelation  function 
of  **x<t>"  when  evaluated  at  =  0  yields  the  average  power 
of  the  process  <E[X2<t)]>.  Thus,  Rf(0>  equals 
the  incoherent  intensity  [25]. 

As  noted  above,  the  power-spectral-density  (PSD)  of  the 
signal  incident  at  the  Bobile  receiver  is  a  primary  factor 
to  be  examined  in  this  study.  According  to  the  Hiener- 
Khintchine  theorem  [25],  the  PSD  of  a  wide-sense  stationary 
process  and  its  autocorrelation  function  fora  a  Pourier 
transform  pair.  Hence, 

S<f>  =  /[Ruu<T)exp<-J2TrfT)  d  t  <8> 


Note  that  the  PSD  and  the  autocorrelation  functions  are 
real  even  functions  of  ”f"  and  ”r”  respectively. 

Furthermore,  the  PSD  nay  be  changed  fro*  a  frequency 
spectrun  to  angular  spectrua  via: 

f  =  <Vo/i)  *  cosa  (9) 

This  latter  distinction  uill  assist  in  interpreting  the 
experimental  results  to  deteraine  the  angular  contribution 
of  the  incoherent  field  relative  to  the  receiver. 

POSER-SPECTRAL  THEORY  OP  PROPAGATION 

The  PSD  theory  of  propagtion  as  presented  here  was  developed 
in  the  early  1970s  by  Bell  Laboratories  [103  -  [123.  It  was 
initially  accepted  as  an  accurate  aodel  based  on  the 
experimental  results  of  Cox  [143-[163.  Further  efforts  to 
confirm  the  theory  have  been  hampered  by  the  very  narrow 
fractional  bandwidth  <2Vo/c;  c  -  speed  of  light)  of  the 
Doppler-spread  signal.  Utilizing  a  novel  approach 
(described  below),  Campbell  [183  was  able  to  obtain  milli- 
hertz  resolution  of  the  received  signal  and  thus  demonstrate 
the  validity  of  the  theory  and  refine  it  for  the  urban 
environment  (some  of  his  results  are  presented  here). 

In  the  Bell  Labs  model,  the  field  at  the  receiver  is  modeled 
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as  the  summation  of  a  number  of  horizontally  traveling  plane 
waves  with  random  amplitudes,  arrival  angles  and  phases  (see 
fig.  3).  The  phases  are  uniformly  distributed  between  0  and 
2*.  The  theory  initially  assumes  that  the  arrival  angles 
are  similarly  distributed.  Given  a  vertically  polarized 
signal,  the  average  (incoherent)  field  is  characterized: 

Ez  =  Eo  ZCn  *  c°s(uct  +  Hnt  +  <|>n)  (10) 

where:  uc  _  transmitted  (radial)  frequency 

“n  ~  doppler  shift  (radial) 

-  random  phase 

Cn  -  normalization  constant  such  that 
<  Cn2>  -  1 

For  large  values  of  N,  the  Central  Limit  Theorem  leads  to 
the  presumption  that  the  field  is  a  narrow-band  Gaussian 
random  process  and,  hence,  a  wide-sense  stationary  process 
C263.  Poliowing  traditional  communications  theory  [27], 
the  field  may  be  further  expressed  in  terms  of  its  in-phase 
and  quadrature  components: 


incoming  wave 


?c<t)  and  Ts(t>  are  therefore  Gaussian 

random  processes  and  possess  zero  mean  and  equal  variance: 


2  m  2 
<TC  >  =  <Tg  > 


<IEZI  > 


<  12) 


where  Tc  and  Ts  signify  the  random  variables  corresponding 
to  the  quadrature  components  at  fixed  "t”  and  are  assumed 
to  be  uncorrelated.  Note  also  that  the  average  power  is 
equal  to  the  variance.  The  envelope  of  Ez<r): 


2  2  .  % 
r  -  <Tc  +  Ta  > 


(13) 


is  Rayleigh  distributed. 

Having  established  the  field  expressions  and  their  relevant 
statistics,  it  now  resains  to  develop  the  predictive 
formulae  for  the  power  spectrum.  Given  a  large  number  of 
scatterers  in  the  the  environsent,  it  is  reasonable  to 
assume  a  continuous  distribution  of  power  incident  on  the 
antenna  within  sose  differential  angle.  Hence,  we  can 
empress  the  differential  received  power  at  the  mobile 
antenna  with  respect  to  arrival  angle  as: 


b  G(a)  p<a)  da 


(14) 


where:  b  -  total  average  power  received  by  isotropic 
antenna. 

G(a)  -  directive  gain  of  receiving  antenna 

p(a)  da  -  fraction  of  total  incoaing  power  within 
da 

Recalling  the  foraula  for  Doppler  shift  as  a  function  of 
arrival  angle  <eq.  1),  and  noting  that  f(a)  =  f(-a),  the 
differential  variation  of  power  as  a  function  of  arrival 
angle  and  frequency  is* 

S(f )  Idf  I  =  b[p(  a)G(  a)  +  p(- a)G<- a)  1 1 d al  (15) 

Since 

a  =  cos  1 - -  (16) 

the  final  expression  for  the  nodel  of  the  PSD  of  the 
received  signal  at  the  aobile  receiver  antenna  (ie.,  the 
situation  as  shown  in  fig.  2)  is: 


14 


S<f>  =  b[p<a)G(a )  +  p<-a)G<-a)]  (17) 

Cf.  -  <f -fc)2ih 

for  If  -  fc  I  <  fm  ,  otherwise  S(f>  =  Q 

This  is  essentially  a  two-dimensional  model,  in  that  all  of 
the  "N"  incident  waves  are  assumed  to  arrive  at  a  near 
horizontal  elevation  angle.  A  more  complicated  model 
accounting  for  a  three-dimensional  distribution  of  arriving 
energy  was  developed  by  Aulin  [28].  However,  the  Campbell 
study  [18],  as  well  as  the  results  of  the  current  research, 
indicate  that  the  simpler  model  presented  above  is  an 
adequate  representation.  Furthermore,  measurements  by  Lee 
and  Brandt  [29]  show  that  the  angle  of  arrival  of  the 
incoming  energy  for  a  mobile  receiver  in  an  urban 
environment  is  concentrated  within  16  degrees  of  the 
horizon. 

e 

The  physical  model  of  eq.  (17)  accounts  for  the  PSD  of  the 
diffuse  component  of  the  received  energy  but  does  not 
include  the  direct  wave.  The  direct  wave  may  be  accounted 
for  in  the  model  by  adding  to  eq.  (17)  the  expression  below. 

B  G<ac  )  fa  -  fc  -  ^  )  +  S  >  (18) 

where:  ac  -  arrival  angle  of  direct  wave 
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fo  *  **  *  cos  ac 

B  -  magnitude  weighting  factor 

Sd(f)  -  PSD  of  amplitude/phase  fluctuations  of 
the  direct  wave 

Equations  <17)  and  (18)  taken  together  represent  the  model 
against  which  the  experimental  results  obtained  in  this 
study  will  be  compared.  Note  that  the  addition  of  the 
coherent  contribution  results  in  a  Rician  probability 
function  for  the  signal  envelope.  A  complete  discussion  of 
the  model  to  include  consideration  of  cross-correlations, 
statistical  moments,  level-crossing  rates,  fading 
properties,  and  coherence  bandwidths,  is  presented  in 
references  1101,  [11],  and  [131. 

The  model  may  be  represented  as  a  frequency  spectrum  as 
above,  or  as  an  angular  spectrum  via  eq.  16.  This  spectrum 
shows  the  distribution  of  intensity  as  a  function  of  arrival 
angle  (fig.  2).  As  the  nature  of  the  scattering  environment 
will  manifest  itself  in  the  expression  for  p(a ) ,  the  angular 
spectrum  presentation  of  the  measured  intensity  gives 
greater  insight  into  the  propagation  mechanism.  Pig.s  4a, b 
show  the  theoretical  normalized  plots  of  the  PSD  as  angular 
spectrum  for  the  case  where  there  is  a  uniform  distribution 
of  scattered  energy  incident  on  the  receiver  antenna  (ie. , 
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p<«)  *  l/(27i')  for  both  a  verticaliy-polari2ed  and  a 
hori2ontally-polari2ed  receiving  half-wave  dipole  antenna 
(transmitter  is  copolari2ed) .  Campbell  [18],  in  refining  the 
model  for  the  suburban  environment,  suggested  a  different 
p(a )  based  on  the  delay-spread  measurements  of  Cox  [15] 
which  showed  that  scatterers  nearest  to  the  moving  receiver 
made  the  greatest  contribution  to  the  received  energy.  Given 
the  scattering  geometry  of  a  vehicle  moving  down  a 
residential  street  with  buildings  on  each  side,  and  noting 
that  the  power  arriving  from  any  one  scatterer  (building) 
should  be  inversely  proportional  to  the  square  of  its 
distance  from  the  receiver,  Campbell  [18]  suggests  that  p(a> 
*  sin2  (x  /it  for  the  urban  environment.  This  results  in  the 
PSD  of  fig.  4c. 

EXPERIMENTAL  APPARATUS 

In  order  to  acquire  meaningful  data  to  verify  the  above 
theory,  the  two  most  significant  problems  that  had  to  be 
overcome  in  the  experimental  process  were  frequency  drift  in 
the  receiver  and  transmitter  and  the  lack  of  an  available 
spectrum  analy2er  with  fractional-hertz  resolution  at  UHP 
radio  frequencies.  The  technique  used  in  this  study  was 
developed  by  Campbell  [18]  and  used  successfully  in  his 
propagation  studies  in  the  suburban  environment.  A  block 
diagram  of  the  equipment  used  is  in  fig.  5.  The  very  stable 


THEORETICAL  ANGULAR  SPECTRUM 

URBAN  STUDr 


sn  3  bids  a 


* 

> 

k 


s 


432  MHZ 


20 


CH  transmitter  is  set  at  432.2  HHz  and  radiates  into  an 
KLH-420- 470-14  antenna  <48  degree  beamuidth).  The  receiving 
antenna  is  a  half-wave  dipole  with  balun.  The  frequency 
drift  of  the  transmitter  was  found  to  be  less  than  0.5  Hz 
over  the  period  that  the  individual  data  runs  were  being 
conducted.  This  resolution,  although  not  ideal,  is 
sufficient  to  discern  the  pattern  of  the  received  PSD. 

The  next  obstacle  uas  the  measured  receiver  drift  of  >  18 
Hz/hour.  This  uas  compensated  for  in  the  following  manner. 
Pirst  the  output  of  a  very  stable  5  HH2  "Vectron"  oscillator 
uas  divided  via  TTL  logic  down  to  100  kHz.  The  resultant 
pulse  has  a  very  high  harmonic  content  into  the  UHP  band. 
This  frequency  "comb”  is  then  inserted  to  one  port  of  a 
hybrid  combiner  at  the  receiver  input.  The  transmitter  is 
tuned  to  a  frequency  about  1200  Hz  below  one  of  the 
reference  harmonics  near  432  HHz.  The  received  signal  is 
then  injected  into  the  other  port  of  the  combiner.  The  lower 
side-band  receiver  is  tuned  so  that  the  reference  harmonic 
at  the  audio  output  is  at  about  1800  Hz  and  the  received 
signal  is  centered  at  600  Hz.  While  the  receiver  will  cause 
both  the  reference  and  the  received  signal  to  vary  over  the 
course  of  a  data  run,  the  difference  frequency  between  the 
two  will  remain  constant  at  1200  Hz.  By  inserting  the  output 
of  the  receiver  (either  directly  or  after  recording 
on  magnetic  tape)  into  the  8th-order  Butterworth  filters 


shown  in  fig.  6  to  separate  the  reference  and  the  signal, 
and  by  applying  then  respectively  to  the  saturating  and 
non-saturating  inputs  of  an  XR-2208  analog  nultiplier  1C, 
the  resultant  output  is  the  PSI3  of  the  received  signal- 
down-converted  and  centered  on  the  audio-range  difference 
frequency  of  about  1200  Hz.  Bhile  any  transmitter  drift 
will  still  be  present  in  the  final  signal,  the  receiver 
drift  will  be  renoved  by  this  mixing  process.  The  specific 
technique  used  here  was  to  play  the  output  of  the  receiver 
directly  into  a  high-quality  "SONY"  TC-D5H  portable 
tape-recorder  for  recording  on  "HAXELL  XL-II"  cassettes 

during  the  data  runs  for  later  processing  in  the  lab. 

* 

The  actual  data  gathering  was  performed  in  the  following 
manner.  Having  selected  a  suitable  environment  in  which  to 
conduct  the  experiment  (see  below),  the  receiver  apparatus 
was  mounted  in  the  back-seat  of  a  1984  Toyota  Corolla  (fig. 
7a).  A  wooden  roof  rack  was  constructed  on  which  to  mount 
the  half-wave  dipole  about  76.2  cm  above  the  roof  of  the 
car  (in  order  to  mitigate  any  effects  of  the  auto  on  the 
received  pattern  and  to  approach  having  the  antenna  in  as 
close  to  a  free-space  environment  as  possible).  See  fig.  7b 
The  car  was  then  driven  to  the  north  end  of  the  selected 
path.  Next,  the  reference  and  signal  frequencies  were 
aurally  cofapa red  with  a  tuning  fork  to  ensure  they  were  at 
about  1800  and  600  Hz.  An  identifying  signal  was  then 


placed  on  the  tape  by  keying  the  reference  on  and  off  with  a 
toggle  switch.  Following  a  one-minute  wait,  the  car  was 
accelerated  up  to  a  steady  84  kph  <+/-  0.8  kph  over  the 
entire  run).  The  reference  was  keyed  off /on  when  the  car 
passed  the  designated  start  point  of  the  run.  The  actual 
length  of  the  path  over  which  constant  velocity  data  was 
taken  was  1.0  km  < vertical/cross-polari2ed  runs)  or  1.6 
km  (horizontal  runs).  The  reference  was  again  keyed 
off/on  at  the  end  of  the  test  path  (the  start  and  stop 
points  of  the  run  were  marked  by  easily  identifiable  terrain 
features) .  The  car  was  then  slowed  to  a  stop  and  turned 
around.  The  complete  procedure  was  repeated  three  more 
tines. 

The  entire  experiment  was  conducted  once  with  both  antennas 
vertically  polarized  (path  length  1.0  kilometers),  once 
horizontally  polarized  (1.6  km  long  path),  and  once  with 
the  transmitter  vertically  polarized  and  the  receiver 
horizontally  polarized  (same  path  as  vertical).  The  data  was 
initially  checked  for  usefulness  by  playing  the  tape  into 
the  f ilters/multiplier  board  whose  output  was  in  turn 
attached  to  a  HP  3682A  Spectrum  Analyzer.  If  accepted,  the 
data  was  played  into  the  HP  5451C  Pourier  Analyzer  for 
storage  and  analysis. 


EXPERIMENT  ENVIRONMENT 


A  nap  of  the  experiment  geography  is  shown  in  fig.  8.  The 
area  is  located  on  the  Keweenaw  peninsula  in  northern 
Michigan.  Data  uas  taken  during  wind-less  days  in  late  May 
and  early  June  when  the  trees  were  in  full-leaf.  The 
transmitter  antenna  uas  located  on  top  of  a  residential  hone 
about  70  meters  from  the  western  shore  of  Portage  Lake.  The 
mobile  receiver  covered  the  routes  shown  in  fig.  8.  The 
road  uas  located  about  8.8  km  from  the  transmitter.  An 
exaggerated  profile  chart  of  the  radio  link  at  about  the 
mid-point  of  the  receiver  path  is  shown  in  fig.  9.  The  total 
elevation  change  over  the  length  of  the  data  runs  uas  about 
5  meters.  At  the  receiver,  the  angle  of  incidence  of  the 
direct  wave  is  less  than  0.4  degrees.  Based  on  the  profile 
and  the  average  tree  height  (see  below),  the  propagating 
radio  wave  first  encounters  the  forest  about  1.6  km  from  the 
receiver.  The  amount  of  forest  penetrated  by  the  direct 
wave  varies  over  the  route  of  the  mobile  receiver  from  just 
less  than  1.0  km  to  about  1.6  km  (as  evidenced  by  the  open 
areas  indicated  on  the  map)  and  averages  about  1.2  km  for 
the  vertical/co-polarized  data  route  and  1.4  km  for  the 
horizontal  route. 

The  forest  uas  surveyed  using  the  random-tree-pairs  sampling 
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■ethod  [30]  -  [32]  (see  Appendix  B  for  a  brief  description 
of  this  technique).  This  method  lends  itself  well  to 
performing  rapid  surveys  of  forested  stands  and  had  been 
utilized  by  the  author  in  the  past.  The  results  of  the 
survey  are  presented  in  table  1.  In  brief,  the  forest 
consists  predominately  of  maple  trees  about  0.7  meters  in 
circumference  and  1.3  meters  tall.  The  coefficient  of 
variation  for  the  tree  diameter  suggests  that  this  is  a 
medium-aged  stand  of  60-70  year  old  trees.  There  are  a  large* 
number  of  "reproduction  class"  trees  (3-10  cm  in  diam.) 
throughout  the  area.  A  light-to-moderate  undergrowth  of 
ferns  and  grasses  is  prevalent  in  the  flatter  eastern 
portion  of  the  area-of-interest  indicated  on  .the  map,  while 
a  much  thicker  and  taller  undergrowth  predominates  over  the 
more-sloping  western  side.  The  crowns  of  the  deciduous 

trees  extend  about  15-30%  down  the  height  of  the  tree  and 
are  characterized  by  an  upright  uide-mouthed-cone  shape  on 
their  undersides  topped  by  a  hemispherical  upper  region.  The 
hemlocks  take  on  a  "christmas-tree”  shape  over  75-80%  of 
their  height. 

The  road  and  the  accompanying  shoulders  and  drainage  ditches 
occupy  a  well-defined  25-meter-wide  gap  in  the  forest.  The 
open  areas  present  in  the  transmission  path  within  1.6  km 
of  the  road  are  characterized  by  fern  growth  60-80  cm  high, 
with  a  smattering  of  young  scotch  pine,  as  well  as  a  few 


TABLE  Bl.  Forest  Survey  Results 


Species 

Average 

Diaaeter 

(ca) 

Dias. 

Stnd. 

Dev. 

Average 

Height 

<») 

Hgt . 

Stnd. 

Dev. 

Relative 

Density 

(4) 

Haple 

23.04 

9.31 

13.64 

2.20 

82 

Birch 

24.23 

7.47 

14.33 

3.17 

10 

Healock 

34.37 

14.46 

13.40 

2.60 

4 

Aspen 

29.68 

9.02 

15.00 

0.43 

3 

Others 

— 

— 

— 

— 

<  1 

Bean  Area  Per  Tree 
(All  types) 

15. 

26  square 

aeters 

Trees  per 
(>  10  cm 

Hectare 
dian. ) 

655 

.14  trees 

per  hectare 

Repro.  Class  Trees 
(3-10  cm  diaa.) 


785.46  trees  per  hectare 


aspen  and  birch  seedlings.  The  transistion  zone  between  the 
open  areas  and  the  aain  stand  of  forest  is  doainated  by 
7-9  aeter  tall  scotch  pines  six  or  nore  aeters  apart.  The 
forestation  pattern  indicated  on  the  aap  is  generally 
accurate  within  the  area  surveyed.  Based  on  the  author's 
brief  survey,  the  pattern  of  forestation  along  the 
transaitted  path  as  it  appears  to  the  nobile  receiver  is 
essentially  a  quasi-stationary  random  process.  It  should 
also  be  noted  that  the  transmitter  antenna  must  radiate  its 
energy  through  a  snail  stand  of  about  30-aeter-high  nixed 
hardwood  trees  before  encountering  open  air. 

PROCESSING  AND  RESULTS 

"Probe"  spectra  are  recordings  of  the  received  signal 
while  the  car  was  waiting  notionless  for  one  minute  prior  to 
accelerating  to  speed  to  start  the  run.  These  are  referred 
to  as  "probes"  since  the  received  PSD's  taken  during  the 
diata  runs  are  the  actual  Doppler-spread  spectra  convolved 
with  the  probe.  Representative  probes  for  each  experiment 
configuration  are  shown  in  fig.  10a  -  c.  The  horizontal  and 
vertical  co-polarized  probes  were  taken  about  0.5  kn  north 
of  the  northern  end  of  the  test  routes  (see  fig.  8)  while 
the  cross-polarized  probe  was  recorded  at  a  point  about  0.5 
ka  south  of  the  southern  end  of  the  data  route. 
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The  received  signal  is  clearly  seen  as  the  sain  peak  in  the 
center  of  the  spectrum  almost  30  dB  higher  than  any  of  the 
other  peaks.  The  recorder  speed  fluctuates  slightly  during 
the  recording  and  playback  process  ("wow  and  flutter")'.  The 
results  of  this  fluctuation  are  easily  seen  as  the  sideband 
peaks  in  the  probe  spectrum.  While  a  totally  "clean”  probe 
<no  sideband  peaks)  is  desired,  the  probe  plots  indicate 
that  the  sidebands  should  have  little  effect  on  the  received 
data.  Furthermore,  the  lack  of  a  spectral  spread  around  the 
received  signal  indicates  the  static  nature  of  the 
environment.  Transmitter  drift  can  also  be  ascertained  by 
examining  the  probes  just  prior  and  following  a  data  run. 

As  noted  above,  the  data  was  processed  using  an  HP  5451C 
Fourier  Analyzer  (see  Appendix  C  for  a  brief  discussion  of 
the  procedures  used).  The  narrow  bandwidth  spectra  were 
obtained  by  using  the  "Band-Selectable  Pourier  Analysis” 
(BSFA)  software  package  that  was  supplied  with  the  machine 
[331  (Bendat  and  Piersol  [34]  discuss  zoom-transform 
techniques  applicable  to  BSFA).  The  resultant  spectral 
resolution  is  0.4  Hz  per  point  (channel).  The  incoming  time 
record  was  multiplied  by  a  single  Hanning  window,  so  the 
equivalent  noise  bandwidth  of  each  channel  is  0.6  Hz  [35]. 

As  noted  above,  frequency  drift  is  on  the  order  of  less  than 
0.5  Hz  so  frequency  uncertainty  within  the  plots  is  bounded 
by  +/-  1  channel.  The  normalized  standard  error  estimate  for 


the  amplitudes  presented  in  the  spectral  plots  is  a  function 
of  the  nusber  of  averages  and  effective  bandwidth  used  in 
saapling  the  data  [343.  During  the  processing  of  the  data 
froa  each  run,  the  BSPA  paraaeters  were  adjusted  to  insure 
that  the  entire  time  record  of  the  run  would  be  averaged 
over  to  produce  the  resultant  power  spectra.  This  resulted 
in  16  averages  for  the  vertically  polarized  and 
cross-polarized  runs  and  22  averages  for  the  horizontally 
polarized  runs  and  the  probe  spectra.  This  yields 
noraalized  standard  errors  for  the  cross-polarized  and 
vertical  spectruas  of  0.354  and  0.302  for  the  probes  and 
horizontal  spectruas.  The  sanpling  distribution  is 
described  by  a  chi-square  distribution  with  32  (vertical  and 
cross  polarized)  or  44  (horizontal  and  probes)  degrees  of 
freedoB  [343.  This  results  in  a  904  confidence  interval  for 
the  true  PSD  of  [0.68  *  MPSD,  1.6  *  HPSD3  for  the  vertical 
and  cross-polarized  results  ("WPSD"  denotes  the  aeasured 
PSD)  and  an  interval  [0.71  *  BPSD,  1.47  *  HPSD3  for  the 
horizontal  results. 

The  received  and  processed  power-spectral-densities  for  each 
of  the  three  polarization  situations  are  presented  in  fig. 

11  -  13.  They  represent  the  single  data  run  out  of  the  four 
for  each  polarization  over  which  vehicle  speed  was  the  nost 
constant  and  the  least  transmitter  drift  occurred.  The 
Doppler-spread  of  the  diffuse  component  is  easily  seen  in 
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all  three  plots  (about  67.2  Hz  wide  yielding  a  vehicle  speed 
of  84.50  kph) .  The  coherent  intensity  peak  is  also  readily 
discernible  at  the  center  frequency  in  each  case. 

The  f ree-space-loss  (FSL)  for  the  line-of-sight  signal  path 
in  the  experiment  is  104  dB  at  432  HHz.  The  dB  scales  of  the 
plots,  which  show  the  excess  path  loss  beyond  PSL,  were 
calibrated  by  carefully  measuring  the  noise  floor  of  the 
receiver,  taking  into  account  the  Hanning  window,  and 
labeling  the  plots  accordingly.  The  noise  floor 
measurements ,  and  hence  the  labeled  plots,  are  accurate  to 
+/-  1  dB.  The  excess  path  loss  is  due  to  the  effect  of  the 
propagating  medium  on  the  radio  wave.  The  -general  down-ward 
slope  of  the  PSD  curve  noticeable  in  the  higher  frequencies 
of  the  horizontal  data  plots  is  due  to  the  tuning  of  the  the 
received  signal  too  close  to  the  upper  edge  of  the  SSB 
receiver's  narrow  pass-band. 

Pigs.  14  -  16  show  only  the  Doppler-spread  bandwidth  portion 
of  the  previous  figures  with  the  horizontal  axis  relabeled 
using  equation  16  to  present  the  data  as  an  angular 
spectrua.  Superimposed  on  the  plots  for  the  co-polarized 
runs  is  the  theoretical  incoherent  PSD  calculated  from  eq. 

17  <ie.,  without  the  direct  wave  contribution)  using  the 

directive  gain  of  the  dipole  for  each  orientation  [361. 
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Pigs.  17  and  18  show  the  spatial  autocorrelations  of  the 
incoherent  intensity,  in  terns  of  nunber  of  wavelengths, 
for  the  co-polarized  runs  (see  Appendix  A  for  derivation 
of  plots).  The  horizontal  scale  was  deternined  via  eq.  7 
and  the  vertical  scale  denotes  a  nornalized  autocorrelation. 
The  snooth  curve  indicated  on  the  plots  is  the  theoretical 
nornalized  autocorrelation  plot  based  on  the  directive  gain 
of  the  two  antenna  polarizations  and  the  equation  for  the 
theoretical  PSD  <eq.  17)  using  a  uniforn  p(a)  for  incoherent 
intensity  arrival  angle. 

DISCUSSION  AND  COMPARISION  HITH  PREVIOUS  HORK 

Pigures  14  and  15  show  the  conparision  of  the-neasured  PSD 
with  the  theoretical  PSD  of  a  uniforn  p(a).  Such  a  p(a) 
indicates  that  there  is  a  uniforn  distribution  of  scattered 
energy  incident  on  the  receiving  antenna  over  a  360  degree 
radius.  Purthernore,  this  inplies  that  the  scattered  wave 
anplitude  is  uncorrelated  with  the  angle  of  arrival.  If  the 
forest  can  indeed  be  nodeled  using  the  uniforn  p(a),  then  a 
long  transnission  path  can  be  represented  at  any  point  along 
the  direction  of  travel  of  the  nobile  receiver  as  a  uniforn 
circle  of  equivalent  scatterers  centered  on  the  antenna, 
despite  the  actual  geonetry  of  the  environnent. 

In  other  words,  scatterers  that  are  at  sone  depth  in  the 
forest  stand  are  just  as  inportant  in  terns  of  received 
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power  as  those  adjacent  to  the  road.  Pig.  19  [18]  shows 
the  contrasting  situation  of  a  non-unifora  p(a>  measured  by 
Caapbell  in  suburban  Seattle,  Hash.  The  fundamental  differ¬ 
ence  between  the  two  environments  which  nay  account  for  the 
distinct  p(a>"s  is  the  more  "fine-grained”  nature  of  the 
forest:  the  trees  are  smaller  and  more  numerous  than  are  the 
buildings.  This  probably  gives  rise  to  higher-order  multi¬ 
ple  scattering  effects  [24]  in  foliage  with  the  scattered 
intensity  being  uniformly  diffused  over  all  arrival  angles 
as  the  wave  propagates  deeper  into  the  forest.  Hence,  one 
must  consider  the  effect  on  total  intensity  of  trees  deeper 
in  the  woods  and  not  just  those  closest  to  the  antenna. 

Lang,  et  al.  [5],  have  postulated  that  as  the  transmitted 
energy  penetrates  deeper  into  the  forest,  the  incoherent 
field  becomes  progressively  more  important.  This  is  evident 
from  figs.  11-13.  A  measure  of  the  ratio  of  coherent- 
to-total-intensity  can  be  obtained  from  these  figures  by 
summing  the  values  of  the  intensity  within  a  feu  channels  of 
the  coherent  peak  and  comparing  it  to  the  total  intensity 
across  the  Ooppler-spread  bandwidth.  Doing  so  yields  a 
value  of  about  0.05  for  the  vertically  co-polarized  case, 
0.42  for  the  horizontal  case,  and  0.21  for  the 
cross-polarized  situation  (values  accurate  to  about  +/-  25% 
of  the  stated  value  for  each  case).  Data  taken  by  Violette, 
et  al,  [6],  indicates  that  for  short  transmission 
paths  through  a  small  number  of  trees,  scatterers  closest  to 


the  receiver  nay  dominate  the  received  intensity  as  in  the 
urban  case  (fewer,  more  discreet  scatterers).  The 
implication  of  the  foregoing  discussion  is  that  the  node  of 
propagation  of  the  transmitted  energy  changes  as  the  depth 
of  penetration  of  the  radiouave  into  the  forest  or  foliage 
is  increased.  Indeed,  the  data  suggest,  as  does  Lang,  et  al 
C51,  that  as  this  depth  increases,  the  coherent  component  of 
the  intensity  may  be  completely  obscured  (especially  for  the 
vertical  co-polari2ed  signal). 


An  important  result  of  any  propagation  study  is  the  observed 
basic  transmission  loss  which  is  usually  expressed  in  terms 
of  an  exponential  "specific  attenuation  rate",  ”a”  (in  dB/m) 
(not  to  be  confused  with  angle  of  arrival  as  used  earlier). 

A  number  of  both  theoretical  and  empirical  studies  have 
attempted  to  examine  and  model  the  rate  for  propagation  in  a 
foliated  environment  (see  references).  Based  on  the 
received  signal  level  of  the  coherent  peats  of  the  PSD  from 
figs.  11  and  12  and  the  average  depth-of-foliage  as  seen  by 
the  receiver,  the  calculated  "a"  for  the  propagation  path 
studied  here  is  0.06  dB/m  for  the  vertically  co-polarized 
case  and  0.04  dB/m  for  the  horizontal  data.  This  is  a 
somewhat  simplistic  representation  of  the  complex 
attenuation  process  over  the  aggregate  length  of  the  signal 
path  and  presumes  that  all  excess  path  loss  is  due  to  the 
foliage  penetration  experienced  within  1.6  km  of  the 
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receiver  (see  fig.  9).  The  possible  adverse  effects  of  the 
coherent  uaves  reflected  from  the  surface  of  the  lake  (fig. 
9),  effects  of  the  slight  elevation  change  over  the  aobile 
receiver  route,  or  froa  other  sources  of  error  are  not 
considered. 

The  specific  attenuation  rates  noted  above  are  an  order  of 
nagnitude  lower  than  many  of  the  aeasureaehts  made  by  others 
in  the  saae  frequency  range  over  shorter  signal  paths  ([71, 
[81,  and  [231  ;  Broun  and  Curry  [31  and  Sass  [91  report 
results  froa  multiple  studies),  but  are  consistent  with  the 
observation  that  the  attenuation  rate  apparently  decreases 
as  depth  of  foliage  penetration  increases  with  the  excess 
path  loss  generally  following  the  pattern  of  fig.  20  [6,91. 
This  figure  is  froa  Sass  [91  (by  peraission  of  the  journal) 
and  shows  excess  path  loss  as  a  function  of  antenna  height 
froa  experimental  tests  of  the  US  Army's  ’’position  locating 
and  reporting  system’’  (PLRS)  in  the  frequency  range  420-450 
BHz.  The  results  from  the  current  study  are  plotted  in  the 
figure  for  coaparision  (polarization  and  forestation 
descriptions  were  not  given  by  Sass).  The  actual  curve  for 
any  specific  situation  will  vary  with  geography,  type  of 
foliage,  antenna  heights,  frequency,  etc. 

The  ratio  of  the  hori2ontal-to-vertical  attenuation  rate  is 
in  near  agreement  with  the  relationship  at  432  MHz 
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FIGURE  20.  EXCESS  PATH  LOSS  CURVE 


hypothesized  by  the  Rayleigh  Effective  Volume  (REV)  model 
for  propagation  in  foliage  developed  by  Broun  and  Curry 
(see  Appendix  D  for  a  brief  description  of  their  nodel) 

[2,33  (see  fig.  D1  for  the  REV  aodel  as  applied  to  a 
hypothetical  forest).  Houever,  the  magnitudes  of  the  - 
individual  rates  are  smaller  than  uould  probably  be 
predicted  by  REV  for  this  study's  environment  as  that  model 
is  only  concerned  uith  direct  uave  propagation  and  does  not 
consider  other  louer-loss  propagation  nodes  that  may  arise 
over  longer  paths  (as  in  fig.  20). 

The  usual  explanation  for  the  phenomenon  described  by  fig. 

20  has  been  [21  -  [5]  the  lateral  uave  propagation  theory 
developed  initially  by  Tamir  and  extended  by  others  [19]  - 
[22].  The  lateral  uave  theory  of  radio  propagation  through 
foliage  models  the  forest  as  a  dissipative  dielectric  slab 
or  half-space,  as  characterized  by  some  effective 
permittivity,  lying  on  an  imperfect  conducting  ground  (the 
earth).  Above  the  slab  lies  free  space  (see  fig.  21).  In 
the  aodel  the  transmitted  uave  leaves  the  source  antenna  and 
some  portion  of  it  is  refracted  at  the  critical  angle  on  the 
air-forest  interface  uhere  it  travels  along  the  boundary  in 
essentially  free-space  (louer-loss)  conditions  continually 
"leaking"  energy  back  into  the  forest,  some  of  uhich 
eventually  reaches  the  receiver.  This  uave  uould  be  much 
less  attenuated  than  the  uave  that  travelled  directly  from 


the  transmitter  through  the  slab  (forest)  to  the  receiver 
and,  hence,  should  dominate  the  received  field  intensity  on 
long  paths.  Tanir  has  shown  [201  that  the  lateral  wave  may 
arise  even  if  the  transmitter  is  not  immersed  in  the  forest 
and/or  more  than  one  stand  of  forest  (separated  by 
free-space)  occurs  in  the  propagation  path;  both  conditions 
which  exist  for  the  current  study. 

Despite  its  explanation  of  the  nature  of  the  excess  path 
loss  curve  of  fig.  20,  recent  studies  reported  by  Violette 
[6]  and  Sass  [9],  which  discuss  measurements  of  parameters 
such  as  impulse  response,  multipath  effects  and  propagation 
delays,  in  addition  to  basic  transmission'loss,  suggest  that 
the  lateral  wave  nay  not  be  a  wholly  accurate  representation 
of  the  propagation  mechanism.  They  propose  instead  a 
"scattering  mode  of  propagation”.  Although  by  no  means 
conclusive,  the  results  of  this  study  also  seem  to  support 
this  contention.  As  noted  above,  one  would  expect  the 
coherent  intensity  to  dominate  the  received  power  spectrum 
given  the  lateral  wave  model.  Also,  if  a  lateral  wave  were 
present,  then  it  seems  reasonable,  based  on  Violette's  [6] 
observation  regarding  scattered  intensity  over  short 
transmission  paths  noted  above,  that  one  night  expect  to  see 
a  definite  correlation  of  intensity  with  arrival  angle 
similar  to  the  suburban  study.  This  would  arise  from  the 
relatively  short  path  the  lateral  wave  would  travel  from  the 
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air-forest  boundary  back  down  to  the  receiver,  but  is 
apparently  not  the  case  in  the  results  of  this  study. 

In  summary,  then,  the  data  presented  here  lend  credence  to 
the  supposition  that  extended-range  coaaunications  through 
forested  areas  is  due  to  a  "scatter  propagation  aechanisa" 
rather  than  a  lateral  uave.  In  view  of  theoretical 
developaents  such  as  the  REV  model  of  Broun  and  Curry  [2,3] 
and  the  work  of  Lang,  et  al  [51,  which  show  that 
scattering  and  absorption  of  radio  waves  are  a  function  of 
the  saaller  (relative  to  about  one-half  wavelength  or  less) 
coaponents  of  the  forest  bio-Bass,  it  is  conceivable  that 
this  nechanism  could  include  a  significant  amount  of  energy 
being  propagated  by  forward-scatter  taking  'place  in  the 
upper  regions  of  the  forest  canopy.  Bore  experimental 
studies  of  the  power  spectrua  and  the  propagation  parameters 
mentioned  above  are  required  for  the  development  of  a 
rigorous  mathematical  description  of  this  proposed  scatter 
propagation  node  T9]. 


Pigs.  13  and  16  are  in  agreement  with  the  well-known 
observation  that  foliage  causes  a  great  amount  of 
depolarization  of  the  transmitted  signal.  The  pattern  of 
the  cross-polarized  PSD  (receiver  horizontally  polarized)  is 
consistent  with  a  uniform  probability  of  energy  arrival 
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angle.  The  specific  attenuation  rate  of  0.054  dB/m  is 


slightly  smaller  than  for  the  vertical-vertical  case.  This 


is  expected  from  the  results  of  the  horizontally 


co-polarized  data,  the  REV  model  of  Broun  and  Curry  [2,3], 


and  other  reports  in  the  literature  [5,6,19-21],  uhich 


document  significantly  smaller  attenuation  in  foliage  of 


horizontally-polarized  uaves  up  to  0.5  -  1  GHz  (above  this 


frequency  range,  no  appreciable  difference  is  noticed  or 


predicted).  This  is  due  to  the  vertical  orientation  and 


dielectric  properties  of  the  tree  trunks  [3,5]  (see  App.  D) 


Purther  analysis  of  the  cross-polarized  data  is  hampered  by 


the  unknoun  effect  of  the  small  patch  of  trees  at  the 


transmitter  site. 


The  spatial  autocorrelation  plots  for  the  tuo-copolarized 


runs  are  presented  in  figs.  17  and  18  along  uith  the  Bessel 


function  curves  of  the  theory  (see  Appendix  A).  The  pattern 


of  the  measured  results  show  a  general  agreement  uith  the 


theory  given  a  uniform  p («).  The  vertical  case  appears  to 


shou  a  shift  from  the  theoretical  plot  by  a  factor  of  0.10 


uavelengths.  This  suggests  that  the  true  receiver  antenna 


pattern  (as  mounted  on  the  vehicle)  is  not  perfectly 


isotropic  in  the  plane  of  reception  (as  the  theoretical  plot 


assumes).  Deviation  of  the  antenna  pattern  from  an 


isotropic  form  uill  stretch  the  theoretical  correlation 


curve  along  the  horizontal  axis  and  uill  dampen  the 
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oscillations  (the  greater  the  deviation,  the  more  pronounced 
the  effects  -  see  the  plot  for  horizontal  case).  Additional 
factors  that  nay  help  explain  the  deviation  from  the 
theoretical  case  for  both  polarizations  are:  1)  the 
environment  is  not  perfectly  stationary,  and  2)  the  p(a)  is 
perhaps  somewhat  more  complicated  than  the  simple  uniform 
expression  discussed  above  (although  the  measurements  do 
suggest  that  this  is  an  adequate  ’’simplification"  for  the 
purposes  of  modeling  the  propagation  channel). 

Lee  [13]  indicates  an  autocorrelation  coefficient  of  0.45 
or  less  is  required  for  the  received  fields  to  be 
uncorrelated.  The  measurements  are  in  accordance  with  his 
observation  that  mobile  radio  signals  are  usually 
uncorrelated  within  a  half  of  a  wavelength.  This  has 
obvious  advantages  in  the  development  of  diversity/combining 
techniques  to  improve  reception  at  the  mobile  receiver 
[10,131. 

CONCLUSIONS 

The  research  presented  here  has  been  a  successful  attempt  to 
verify  the  stochastic  power-spectral-density  theory  of  radio 
propagation  as  applied  to  a  forested  environment.  Using  a 
novel  technique  to  achieve  fractional-hertz  resolution  of 
the  Doppler-spread  received  power  spectrum,  the  resultant 
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data  taken  at  432  WH2  over  long  transmission  paths  through 
the  foliage  shows  reasonable  agreement  with  the  model  as 
developed  by  Bell  Labs,  given  a  uniform  probability 
distribution  for  the  angle  of  arrival  of  incident  energy  on 
the  receiving  antenna.  The  data  also  provide  a  possible 
insight  into  the  mode  of  propagation  prevalent  on  the 
inter-foliage  path.  The  results  tend  to  support  the 
contention  of  Hufford,  et  al  [4].  Violette,  et  al  [6],  and 
Sass  [9]  that  the  transition  in  propagation  in  this  medium 
from  a  direct  wave  to  a  lower-loss  mode  over  longer  path 
lengths  may  be  better  characterized  by  a  forward-scattering 
mechanism  than  by  a  lateral  wave  model.  The  measured  ratios 
of  coherent- to-total  intensities  reported  here,  which  show  a 
large  amount  of  incoherent  (scattered)  intensity  present  in 
the  received  signal,  together  with  the  discovery  of  the 
uniform  p(oO  for  the  intensity  and  the  significant 
depolarization  of  the  signal,  suggest  that  high-order 
multiple-scattering  theories  (see  [24]),  such  as  the 
approach  of  Broun  and  Curry  [2,3],  may  hold  the  key  to 
accurate  mathematical  descriptions  of  the  propagation 
phenomenon  in  this  medium. 

The  results  of  this  research  have  definite  implications  for 
the  use  of  diversity  techniques.  The  rapid  decorrelation  of 
the  signal,  as  evidenced  in  figs.  17  and  18,  suggests  the 
use  of  space  diversity  to  improve  reception  at  mobile 
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receivers.  The  required  antenna  spacing  of  one  -  half  - 
wavelength  (evident  from  the  spatial  correlation  plots)  is 
certainly  obtainable  on  most  automotive  platforms  for 
frequencies  in  the  upper  VHP  and  UHP  bands  (see  Lee  [133  for 
information  on  combining  technologies  necessary  to  implement 
space  diversity).  As  the  receiver  path  will  not  always  be 
oriented  orthogonally  to  the  transmitted  signal  as  in  this 
study,  further  research  will  be  needed  to  determine  the 
optimum  placement  of  the  antenna  group  on  the  vehicle. 

The  significant  depolarization  of  the  transmitted  signal,  as 
evidenced  by  figs.  13  and  16,  suggests  that  the  use  of 
polarization  multiplexing  is  not  a  viable  means  by  which 
to  achieve  efficient  use  of  limited  frequency  resources  in 
applications  requiring  transmission  through  foliage. 
Purthermore,  the  results  of  this  study  also  confirm  that 
horizontal  polarization  is  the  "po.larization-of-choice"  be¬ 
low  500  MHz  in  foliage  (see  Appx.D).  The  data  gathered  here 
suggests  that  horizontally  polarized  signals  experience  less 
excess  path  loss  and  retain  a  larger  percentage  of  the 
signal  in  the  coherent  component  than  do  vertically 
polarized  signals  over  the  same  foliated  transmission  path. 

It  is  recommended  that  the  specific  experiment  described 
above  be  repeated  over  the  same  data  paths  at  other  times 
during  the  year  to  assess  the  effects  of  seasonal  changes  in 


foliage.  While  the  geography  uill  not  allou  data  to  be 
compared  from  a  multitude  of  signal  paths  at  a  variety  of 
forest  depths,  the  existence  of  an  unimproved  road  just  to 
the  west  of  the  road  utilized  in  this  study  uill  provide  a 
signal  path  that  has  about  one-half  of  the  foliage 
penetration  reported  here.  It  may  also  be  instructive  to 
gather  data  over  both  of  these  roads  utilizing  various 
antenna  heights. 

The  specific  techniques  described  in  this  thesis  can 
certainly  be  utilized  to  study  propagation  in  other  locales 
and  in  different  mediums.  Purther  research  at  multiple 
frequencies  of  interest  and  over  a  number  of  paths  varied  i 
both  type  of  foliage  and  depth  of  penetration  are  required 
to  provide  a  definitive  description  and  facilitate 
development  of  theoretical  predictive  models  for  UHP 
propagation  in  forests. 
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APPENDIX  A:  DERIVATION  OP  THEORETICAL  PLOTS 

POSER  SPECTRAL  DENSITIES 

The  PSD  plots  of  fig.4a-e  were  obtained  in  the  following 
■anner.  First,  given  a  uniform  p<a>  and  an  isotropic  - 
antenna  (as  in  the  vertical  co-polarized  case),  the 
expression  for  the  normalized  PSD  in  eq.  17  simplifies  to: 

f"fc  2  J* 

s<f)  =  l/ti  -  < — (ad 

u 

This  yields  fig.  4a  when  plotted  as  an  angular  spectrum 
(using  eq.  16)  and  presented  in  dB. 

The  power  pattern  in  the  elevation  plane  of  a  half-wave 
dipole  is  [363: 

Ceos 2< ---case  >3  /  sin2<0)  (A2) 

2 

For  the  antenna  orientation  in  this  experiment,  a  =()  (see 
fig.  7b).  The  plot  of  eq.  A1  multiplied  by  eq.  A2  results 
in  fig.  4b,  the  plot  for  the  horizontal  polariztion. 

In  the  suburban  study,  Campbell  C183  used  a  small  dipole 
with  the  following  pattern: 


A2 

2  , 

G<a)  =  sin  <a>  <A3> 

This  results  in  the  expression  for  the  nornalized  PSD: 

S<f>  =  Cl  -  <- - —  <A4) 

and  plotted  as  an  angular  spectrun  in  fig.  4c. 


THEORETICAL  SPATIAL  CORRELATIONS 


The  theoretical  spatial  correlations  (of  incoherent 
intensity  only)  were  obtained  as  follows.  Since  the 
aeasured  PSD  is  a  real  even  function,  its  autocorrelation 
function  is  expressed  [10]: 

Rn)1< t)  *  2  /WS<f)  cos [2  tt ( f  -  f  >tldf  (A5) 

UU  4-1J,  C 

Substituting  in  eq.  A5  the  expression  for  the  PSD  for  the 
vertical  case  <eq.  Al),  and  noting  froB  eq.  16  that: 

df  =  -  fm  sin  <«>  da  (A6) 

f'f„  2  V 

sin(a )  =  [1  -  ( - 


And  sisplifying  the  resultant  equations  while  changing  the 
Halts  of  integration  via  eq.  16,  results  in  an  expression 


for  the  normalized  autocorrelation  of  the  vertical 
copolarized  case: 

2  * 

R  <t)  =  -  J P  cos[27rf  (cos  a)t]  da  <A7) 

UU  7T  •  iri 

which  is  also  a  Bessel  function  of  the  first  kind  with  order 
0.  To  obtain  the  plot  in  fig.  17,  eq.  A7  was  evaluated  at 
various  values  of  ”t"  utilizing  a  ’’BASIC”  language  computer 
program  incorporating  the  ’’Simpson's  Rule”  algorithm  of 
numerical  analysis  [373.  The  results  were  accurate  to  at 
least  4  decimal  places  (as  compared  to  published  tables). 

The  theoretical  plot  for  the  horizontal  case  (fig.  18)  was 
done  by  the  same  program,  however  the  integrand  was  eq.  A7 
multiplied  by  eq.  A2  (antenna  pattern  with  0  =  a)-  The 
temporal  autocorrelation  results  obtained  by  the  above 
process  could  be  converted  to  spatial  correlations  via  eq. 

7  with  the  distance  expressed  in  wavelengths. 

MEASURED  SPATIAL  CORRELATIONS 

In  order  to  directly  compare  the  measured  data  to  the 
theoretical  incoherent  spatial  correlations,  it  was 
necessary  to  remove  the  coherent  component  of  the  PSD  and 
then  take  the  inverse  Pourier  Transform  to  obtain  a 
"processed”  spatial  correlation.  Using  the  point-by-point 
channel  fill  capability  of  the  HP  5451C,  the  few  channels  of 


the  coherent  component  portion  of  the  168-channel-uide 
Doppler-spread  spectrun  <the  "peaks"  of  the  PSD  plots)  were 
replaced  in  the  data  block  uith  values  siailar  to  the 
adjacent  non-coherent  channels.  The  resultant  data  block 
uas  coapared  both  visually  and  nuaerically  with  the  original 
to  insure  a  "proper"  fit  of  the  "aassaged"  data,  ie. ,  a 
data-block  that  uould  have  occured  if  no  coherent  peak  uas 
present.  This  "processing"  technique  uas  used  for  each  of 
the  three  polari2ation  situations  exaained  in  the  study. 

The  HP  5451C  can  only  perfora  an  inverse  Pourier  transfora 
on  the  full-uidth  512  channel  data  block  of  figs.  11  -  13. 

To  obtain  a  transfora  of  the  168-channel  Doppler-spread 
portion  of  the  "processed"  data-block,  it  uas  necessary  to 
print  out  the  channel  values  and  transfer  then  by  hand  to  a 
personal  coaputer  disk  file  for  further  calculation.  A 
"BASIC"  language  prograa  uas  uritten  using  the  follouing 
equation  to  perfora  the  inverse  transfora  [33],  [34]: 

1*7 

f<n  At)  =  £  P<a  Af)  cos  (2  Trna/168)  (A8) 

rna  o 

This  equation  treats  the  data  as  having  a  data-block  size  of 
168  uhile  aaintaining  the  0.4  Hz  spacing.  Since  the 
relationships  betueen  block-size  N,  Af,  &t,  and  T  (total 
length  of  the  tiae  record)  are  as  follous  [34]: 


T  =  1/  Af 


<A9) 


WBHUIM1  JWWfWIUifWMI  WHM  (WffiflVIlHIWWn  W* 


A5 

;  At  =  T/N  <A9) 

the  resultant  points  of  the  autocorrelation  function  from 
eq.  A8  are  spaced  at  one-half-wavelength  intervals  and  the 
plots  can  be  labeled  accordingly.  The  autocorrelation  is 
normalized  by  dividing  all  values  by  Ruu(0),  which  is  the 
total  power  in  the  PSD  <ie,  integration  of  the  PSD  curve). 


APPENDIX  B:  RANDOM-PAIRS  METHOD  OP  FOREST  SURVEY 


The  randoa-pairs  aethod  is  a  proven  technique  for  Baking 
fairly  rapid  and  accurate  surveys  of  forest  stands  [30-1  - 
[321.  It  is  based  on  the  preaise  that  the  nuaber  of  trees 
per  unit  area  can  be  calculated  froa  the  average  distance 
between  trees.  If  one  considers  the  trees  in  a  forest  to 
occupy  the  centers  of  regular  contiguous  hexagons,  then  the 
aean  area  per  tree  (HA)  is  : 

HA  =  (0.5  d>2  *  3.4644  <B1> 

or  HA  =  d2  *  0.8661  =  (0.93  d>2 

where  ”d”  is  the  average  distance  between  trees.  Dividing 
HA  into  the  unit  area  of  interest  yields  the  trees  per  unit 
area. 

The  field  conduct  of  the  survey  aethod  is  not  quite  as  easy 
as  it  first  aay  seen.  Siaply  taking  trees  at  randon  in  the 
woods  and  aeasuring  the  distance  to  their  nearest  neighbors 
will  result  in  a  density  figure  that  is  greater  than  the  true 
figure.  This  results  froa  aeasuring  the  distance  between 
the  two  closest  trees,  not  between  trees  that  are  an 
average  distance  apart. 


B2 

To  insure  valid  measurements ,  data  on  the  trees  are  taken  as 
shown  in  fig.  Bl.  Sample  points  are  picked  in  the  uoods  by 
some  method  useful  to  the  situation  as  determined  by  the 
surveyor.  The  technique  used  by  the  author  was  to  take 
points  50  -  60  paces  apart  on  East-Best  compass  lines 
through  the  area  of  interest  in  fig.  8.  A  line  is  then 
drawn  (by  stretching  a  string)  from  the  point  to  the  closest 
tree.  Another  line  perpendicular  to  the  first  is  then 
placed  through  the  sample  point  (as  in  fig.  Bl).  This  line 
then  forms  a  180  degree  "exclusion  angle"  with  its  vertex  r 
at  the  sample  point.  Other  exclusion  angles  such  as  160 
degrees  may  be  used  (so  long  as  calculations  are  adjusted  as 
noted  belou)  but  the  180  degree  angle  is  easiest  to  use  in 
the  field.  It  is  called  the  "exclusion  angle”  because  all 
trees  that  fall  inside  of  it  may  not  be  considered  uhen 
looking  for  the  next  tree  in  the  "random  pair".  This  second 
tree  of  the  pair  is  taken  as  the  tree  that  lies  closest  to 
the  line  forming  the  exclusion  angle  but  outside  of  the 
exclusion  zone. 

Initial  comparision  of  the  random-pairs  technique  using  the 
formula  in  eq.  Bl  to  more  comprehensive  forest  surveying 
procedures  showed  that  the  former  yielded  biased  estimates 
of  HA  [313.  It  was  found  that  the  angle  of  exclusion 
corresponded  to  a  specific  correction  factor  in  the 
calculation  for  the  HA  such  that: 


"True"  BA 


<C.F.  tines  "d">2 


<B2> 


The  correction  factor  (C.P.)  for  various  angles  of  exclusion 
are  as  follows  <180  degree  angle  was  used  in  this  research): 


1GLE 

C.F. 

0 

1.67 

30 

1.43 

100 

1.03 

140 

0.89 

160 

0.83 

180 

0.80 

• 

In  practice,  this  Method  is  used  to  sample  trees  in  the 
forest  with  diameters  at  breast  height  <DBH)  of  4  or  more 
inches.  Trees  snaller  than  this  are  considered  as 
"reproduction  -  class"  specimens  and  were  surveyed  in  this 

report  by  taking  the  average  amount  found  in  a  number  of 
randomly-selected  2600-square-foot  plots  in  the  area-of- 
interest  and  translating  that  to  a  figure  for  average  number 
per  unit  area  <see  table  1). 

Tree  heights  were  estimated  using  a  "HERRITT"  hypsometer. 
This  is  simply  a  trigonometrically-calibrated  stick  held 
vertically  out  from  the  surveyor's  body  at  about  three- 
quarters  of  an  arms-length  while  he  is  standing 
a  prescribed  distance  from  the  tree.  The  tree  height  is 
then  read  by  noting  where  the  top  of  the  tree  registers  on 
the  stick. 


APPENDIX  C:  HP  5451C  FOURIER  ANALYSIS 

The  HEHLBTT-PACKARD  5451C  Pourier  Analyzer  provides 
■ulti-channel  analysis  of  DC  -  25  kH2  tine  signals.  Its 
■ain  components  are  a  12-bit-uord  Analog-Digital  converter, 
a  specialized  60K  Processor,  a  key-board-entry  Control  Unit, 
an  oscilloscope  display,  2.6H  hard-disk  drive,  and  a 
magnetic  tape  drive  unit,  along  uith  graphics  terminal  and 
printers  and  plotters.  It  manipulates  data  uith  16-bit 
uords  that  are  collected  in  data  blocks  of  up  to  4096  uords. 
The  reader  is  referred  to  reference  C331  for  complete 
information  on  capabilities  and  operation. 

The  relationships  betueen  the  time  and  frequency  domain 
parameters  of  the  sampled  data  in  its  digital  form  are  as 
follous: 

Time  domain  sample  At  =  l/2Fmax 

point  increment  (Cl) 

=  T/N 

where:  Fmax  -  highest  frequency  in  time  record 

(in  pouers  of  2) 

T  -  length  of  time  record  (as  set  on  ADC) 


N  -  block  size  (as  set  on  keyboard) 


Frequency  increment 


1/T 


<C2) 


Af  = 


=  2Paax/N 


The  Pourier  Analyzer  uses  the  following  equations  to  perfor 
Pourier  Transfons  [33],  [34]: 

F(a  Af)  =  <1/N)  *Ef  (n  At)  exp<-j2Tr  na/N) 

too 

(C3) 

f(n  At)  =  rP(a  Af)  exp<+]2  una/N) 

n>o 

fihen  taking  an  inverse  transfora  (second  equation),  the 
analyzer  assumes  a  symmetrical  spectrua  and  therefore 
coaputes  a  pure  real  tiae  series  [33]. 

Processing  a  single  tiae  record  (of  length  "T")  results  in 
an  inconsistent  estiaate  of  the  PSD  (noraalized  standard 
error  equals  1;  see  [34]).  A  consistent  estiaate  can  only 
be  achieved  by  averaging  a  nuaber  of  tiae  records  (see 
text).  The  nuaber  of  averages  available  for  processing  the 
data  of  this  study  was  liaited  by  the  total  length  of  the 
individual  data  runs  (see  below). 


Using  the  basic  hardware  configuration  of  the  5451C,  the 
best  frequency  resolution  we  might  obtain  in  this  study  was 


C3 


1  Hz  per  channel  (based  on  Paax  of .  2048) .  However,  the 
Band-Selectable-Pourier-  Analysis  software  present  on  the 
disk  operating  systew  allows  greater  channel  resolution  and 
concentration  of  processing  about  a  selected  band  of 
frequencies  within  in  the  original  sawple  record.  This. is 
obtained  through  a  fora  of  coaplex  aodulation  where-by  the 
saapled  band-liaited  frequency  records  (as  selected  by  the 
specific  frequency  band  the  user  wishes  to  concentrate  the 
processing  on)  are  translated  to  base-band.  This  lowers  the 
Paax  in  the  record  and  allows  increased  frequency  resolution 
[341. 

The  HP  5451C  was  set-up  as  follows  to  process  the  data  for 
this  report.  The  "TOTAL  TIHE"  (length  of  a  saaple  record) 
was  set  to  1  second  on  the  ADC  and  only  channel  A  was  used. 
The  BSPA  prograa  was  called  via  key-board  coaaands.  Analog 
through-put  block  size  was  set  to  4096.  BSPA  block  size  was 
set  to  1024  (saxiauw  available;  yields  512  channels  in  the 
frequency  dowain).  The  number  of  tiwe  dowain  records  that 
could  be  saapled  were  deterained  by  the  following  equation 
<BS  =  blocksize). 

BSPA  BS 

ft  Thru-put  Records  =  Zooa  Power  * -  *  #  averages  (C4) 

(of  length  set  T-put  BS 

on  ADC) 


Zooa  Power  =  Paax/ (band-width  to  be  exaained) 


(C5> 


Jj 


C4 

The  frequency  band  ue  wished  to  concentrate  on  was  the  150 
Hz-wide  band  centered  on  the  ^1200Hz  center  frequency  of  the 
out-put  of  the  drift-correction  board  (see  fig.  5  and  text). 
The  closest  the  processor  could  coae  to  our  desired  150  Hz 
bandwidth  was  204.8  Hz.  This  yielded  a  zoos  power  of  ten, 
and  Halted  the  processing  to  either  16  or  22  averages  in 
order  to  include  data  froa  the  entire  length  of  each  data 
run  (41.5  seconds  total  recording  tine  for  each-  vertical  and 
cross-polarized  run  and  55  seconds  for  each  probe  and 
horizontal  run  -  see  text;  aore  averages  also  decrease  the 
standard  noraal  error).  This  results  in  a  Af  of  0.4  Hz 
for  the  BSPA  processed  data. 


t| 


APPENDIX  D:  RAYLEIGH-EPFECTI VE-VOLUNE  MODEL  [23,  [33 


The  following  is  a  brief  presentation  of  Broun  and  Curry's 
REV  model  of  direct  uave  propagation  in  forests  [23,133. 

The  model  is  based  on  the  observation  that  the  uoody 
components  of  a  forest  comprise  less  thqnO.5%  of  the  total 
volume  (despite  uhat  one  may  intuitively  think  from  just 
looking  at  a  forest;  a  case  of  not  being  able  to  see  the 
unoccupied-air-space  for  the  trees!).  This  allowed  Broun 
and  Curry  to  use  the  Foldy-Twersky  theory  of  uave 
propagation  thru  a  sparse  concentration  of  discreet  random 
objects  [24,  vol.23.  Noting  that  attenuation  and  scattering 
of  radio  waves  is  associated  uith  the  smaller  (relative  to 
one-half  wavelength)  components  of  the  forest,  they 
developed  a  fairly  simple  algorithm  to  determine  the 
effective  fractional  volume  of  wood  ( /> Vp>  necessary  for  the 
equations  of  Poldy-Tuersky  (see  [23). 

The  basic  equations  of  the  REV  model  are  as  follows: 

k  -  ko  (  1  +  (/>Vp)  (  e  -  1) (q)tl-  J0.13)  (Dl) 

Im  (k)  =  -ko  (/>VpXq)[  er"  +  (  e  '  ~  1)/103  (D2) 

where:  k  -  uave  number  in  foliage 


D2 


ko  -  free  space  uavenuaber 


-  complex  permittivity 


(  £r"  -  imaginary;  £r'  -  real) 


-  polarization  factor 


=  0.25  for  vertical  polarized  waves  <  500MHz 
=  0.125  horizontal  polarized  waves  <  SOOIlHz 
=  0.167  any  polarization  for  freq  >  50011Hz 


The  difference  in  the  factors  for  the  polarizations  results 
froa  the  vertical  orientation  of  the  tree  trunks,  which 
dominate  the  constituents  of  <yoVp)  below  500  MHz.  Above 
this  frequency,  most  of  the  trunk  material  drops  out  of  the 
effective  volume.  The  branches  and  leaves'  become  the  major 
components  but  are  much, more  randomly  oriented. 


The  real  part  of  the  dielectric  constant  of  wood  is  centered 
around  a  value  of  about  50  while  the  loss  factor  (imaginary 
part  of  the  complex  permittivity)  varies  as  a  function  of 
frequency  and  season.  Both  are  variable  with  tree  type. 


The  plot  of  fig.  D1  results  from  eq.  02  using  a  value  of  50 
for  the  real  part  of  the  permittivity  and  values  for  the 
loss  factor  as  taken  from  fig.  2  of  reference  [31.  It  is 
unknown  what  the  effective  value  of  the  complex  permittivity 
or  the  effective  volume  1/oMp)  are  for  the  forested 
environment  used  in  this  study. 


Sffvivi**w 


It  should  be  reiterated  that  this  model  does  not  a'  apt  to 
explain  the  lou-loss  propagation  mode  phenomenon  of  fig.  20, 
but  rather  is  concerned  only  with  the  direct  wave  that 
propagates  directly  from  the  transmitter  to  the  receiver 
(see  fig.  21).  Broun  and  Curry  accept  the  lateral  wave 
model  to  explain  fig.  20. 


2/2 


AD-A171  431 
UNCLASSIFIED 


INVESTIGATION  OF  THE  UHF  MOBILE-RADIO  PROPAGATION 
CHANNEL  IN  A  TEMPERATE-ZONE  FORESTED  ENVIRONMENTS) 

ARHV  MILITARY  PERSONNEL  CENTER  ALEXANDRIA  VA  U  WHITE 
13  AUG  8£  F/G  17/2. 1 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1%3-A 


CURRICULUM  VITAE 


Kenneth  Janes  White  was  born  13  September  1956  in  Kalamazoo, 
Michigan.  He  was  graduated  from  Port  Huron  Northern  High 
School,  Port  Huron,  Ml,  in  1974.  He  studied  forest  and 
stream  ecology  and  Military  science  at  The  University  ol 
Michigan,  Ann  Arbor,  MI,  graduating  in  1978  uith  a  Bachelor 
of  Science  degree  in  Natural  Resources  (BSNR).  In  July  of 
that  year,  he  was  commissioned  as  a  Second  Lieutenant  in  the 
US  Army  Signal  Corps.  After  receiving  training  in  combat 
communications  at  PT  Gordon,  GA  and  PT  Sill,  OK,  LT  White 
served  for  five  years  in  the  8th  (US)  Infantry  Division 
(Mechanized)  in  Best  Germany.  His  duty  assignments  included 
Communications  Staff  Officer  for  the  division's  Air  Defense 
Artillery  battalion.  Intelligence  Officer  and  later 
Operations  Officer  (S2/S3)  for  the  8th  Signal  Battalion,  and 
commander  of  the  signal  battalion's  Forward  Communications 
Company.  Upon  his  return  to  the  states,  CPT  Hhite  attended 
the  Dept,  of  Defense's  senior  communications  school,  the 
Telecommunications  Systems  Staff  Officer  Course  at  Keesler 
Air  Porce  Base,  MS.  He  was  subsequently  selected  for 
advanced  civilian  schooling  and  attended  Michigan 
Technological  University  in  1985-86  to  pursue  a  Master's 
degree  in  Electrical  Engineering.  His  follow-on  duty 
assignment  is  with  the  Arny  Space  Program,  Los  Angeles  Field 
Office  ,  LA  Air  Porce  Station.  He  and  his  wife  Lynn  are  the 
proud  parents  of  two  sons,  Leonard  and  Darren. 


